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FOREWORD 

This  report  constitutes  the  Final  Report  performed  under 

contract  Nonr  1866(49)  at  the  Gordon  McKay  Laboratory.  The  research 

carried  out  for  Project  DEFENDER  was  concerned  with  the  Optical 

Properties  of  Materials  at  Very  High  Light  Intensities.  This  work  is 

closely  related  to  research  carried  out  in  the  general  field  of  Nonlinear 

Optics  supported  under  the  continuing  contracts  Nonr  1866(28)  and 

Nonr  1866(16).  Some  of  the  research  reported  here  received  joint  support 

under  the  three  contracts.  This  Final  Report  reviews  the  entire  period  of 

the  contract  Nonr  1866(49),  but  more  specifically  the  period  from 

June  1,  1964,  through  June  1,  1965.  The  following  personnel  has  been 

associated  with  the  project  during  the  final  year: 

Professor  N.  Bloembergen  Project  Director  and  Principal 

Investigator 

Dr.  J.  Ducuing  Research  Fellow 

Dr.  R.  K.  Chang  Graduate  Student 

Mr.  P.  Lallemand  Graduate  Student 

In  the  previous  period,  Dr.  J.  vander  .Ziel,  research  fellow,  and 
Mr.  L.  Malmstrom,  a  graduate  student,  have  also  been  associated  with  the 
project. 
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Since  all  essential  results  obtained  under  this  contract  have 
already  been  reported  in  the  published  literature,  this  Final  Report  consists 
only  of  a  Summary  to  which  the  reprints  of  published  papers  are  attached 
as  appendices.  Related  papers,  which  did  not  receive  support  under  contract 
Nonr  1866(49),  have  not  been  included. 


SUMMARY 


The  two  major  experimental  efforts  partially  supported  by 
this  contract  were  concerned  with  A)  The  second  harmonic  production 
of  light  in  semiconductors  and  B)  The  stimulated  Raman  effect  in  liquids. 

A)  The  complex  nonlinear  susceptibility  describing  the  second 
harmonic  generation  of  light  in  semiconducting  I1I-V  and  II-VI  compounds 
has  been  measured  as  a  function  of  frequency.  Both  the  real  and  imaginary 
parts  have  been  determined.  The  determination  of  the  phase  of  the  complex 
susceptibility  has  been  accomplished  for  the  first  time.  This  novel  method 
is  described  in  Appendix  I.  Another  new  result  is  the  variation  of  the  non¬ 
linear  susceptibility  as  a  function  of  frequency  in  the  same  sample.  The 
experimental  results  are  shown  in  Appendix  II.  Further  experimental 

and  theoretical  investigations  in  this  new  field  of  nbnlinear  spectroscopy 
of  solids  will  be  carried  out  under  other  contracts. 

B)  A  Raman  amplifier  cell  technique  has  been  developed  to 
give  a  more  reliable  determination  of  the  stimulated  Raman  gain  in  fluids. 
This  method  is  described  in  Appendix  III.  The  work  in  Appendix  IV  was 
presented  at  the  International  Conference  on  the  Physics  of  Quantum 
Electronics,  held  in  June  1965,  in  San  Juan  (P.  R. ) .  The  most  striking 
result  is  that  the  gain  per  unit  length  in  a  cell  filled  with  Raman  fluid 
changes  as  the  laser  beam  progresses  into  the  medium. 
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Furthermore,  the  initial  gain  per  unit  length  depends  on  the  distance  of 
separation  between  the  Raman  laser  oscillator  and  the  amplifier  cell. 
Subsequent  work  with  the  same  equipment  has  shown  that  the  intensity 
distribution  in  the  laser  beam  changes  as  it  traverses  through  the  fluid. 
Experiments  with  different  fluids  in  two  Raman  cells  in  series  have 
demonstrated  that  this  effect  is  due  to  self- focusing  of  the  laser  beam. 

There  is  a  transverse  gradient  in  the  index  of  refraction,  because  this  index 
is  dependent  on  the  intensity.  The  sell* focusing  effect  is  entirely  distinct 
and  separable  from  the  stimulated  Raman  effect.  It  is  most  proryounced 
in  fluids  with  anisotropic  molecules  and  is  well  correlated  with  the  aniso¬ 
tropy  in  the  optical  polarizability  of  these  molecules.  The  effect  is 
described  in  Appendix  V.  The  discrepancies  of  one  or  two  orders  of 
magnitude  in  the  experimental  and  theoretical  values  of  the  Raman  gain 
in  certain  fluids  are  thus  largely  explained  by  the  creation  of  filaments  of 
high  intensity  in  the  laser  beam  dje  .to  self -focusing.  These  experiments 
have  obviously  important  consequences  for  the  design  of  Raman  lasers 
and  the  passage  of  high  intensity  laser  beams  throjugh  media  with  aniso¬ 
tropic  molecules.  Further  experimental  work  on  factors  determining  the 
size  and  distribution  of  the  bright  filaments,  on  the  influence  of  focusing 
on  the  stimulated  Brillouin  effect,  on  the  magnitude  of  focusing  effects  and 
Raman  gain  in  fluids  with  isotropic  molecules,  etc. ,  is  obviously  suggested 
and  will  be  carried  out  under  other  contracts. 


* 


APPENDIX  1 


RELATIVE  PHASE  MEASUREMENT  BETWEEN  FUNDAMENTAL  AND  SECOND- HARMONIC  LIGHT* 

R.  K.  Chang 

Gordon  McKay  Laboratory,  Harvard  University,  Cambridge,  Massachusetts 

J.  Ducuing 

Physics  Department,  Massachusetts  Institute  of  Technology,  Cambridge,  Massachusetts 

N.  Bloembeigent 

Departments  of  Electrical  Engineering  and  Physics,  University  of  California,  Berkeley,  California 

(Received  14  June  1965) 

The  nonlinear  susceptibility  describing  the  linear  susceptibility, 

second-harmonic  generation  of  light  is  real,  The  experimental  arrangement  is  shown  in 

when  the  medium  is  transparent  at  both  the  Fig.  1.  A  linear  polarized  laser  beam  from 

fundamental  and  the  harmonic  frequency.  When  a  Q-switched  ruby  laser  enters  an  evacuated 

the  medium  is  absorbing  at  either  or  both  of  box  and  generates  second-harmonic  radiation 

these  frequencies,  this  nonlinear  susceptibility  by  reflection  from  a  nonlinear  mirror  of  a 
is  a  complex  quantity.1  The  nonlinear  polar-  crystal  with  43 m  symmetry,  whose  complex 

ization  has  a  phase  shift  with  respect  to  the  nonlinear  susceptibility  must  be  determined, 

fundamental  field.  Measurement  of  the  second-  The  laser  beam  subsequently  generates  an 

harmonic  intensity  only  determines  the  abso-  additional  second-harmonic  field  in  a  potassium 

lute  value  of  the  nonlinear  susceptibility.’  For  dihydrogen-phosphate  (KDP)  platelet  which 

some  IH-V  and  II-  VI  compounds,  which  are  serves  as  a  reference  signal.  The  second- 

absorbing  at  the  second-harmonic  frequency,  harmonic  fields  generated  in  the  two  samples 

such  measurements  i.ave  been  carried  out  in  will  have  a  definite  phase  relationship  with 

reflection. 2»3  In  this  note  an  experiment  is  de-  respect  to  each  other,  as  each  has  a  specific 

scribed  which  determines  the  phase  of  the  non-  relationship  to  the  phase  of  the  fundamental 
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FIG.  1.  Diagram  of  the  experimental  apparatus  to  determine  the  relative  phase  of  the  second-harmonic  field,  as 
described  in  the  text.  Each  laser  pulse  is  monitored  by  the  harmonic  production  in  quartz. 


field  squared.  The  total  second-harmonic  out¬ 
put  is  the  result  of  the  interference  between 
the  two  second-harmonic  beams.  This  inter¬ 
ference  can  be  made  visible  by  admitting  dry 
air  into  the  box.  The  dispersion  in  air  changes 
the  relative  phase  between  the  fundamental 
and  second  harmonic  in  the  5. 5- cm  path  between 
the  reflection  spot  on  the  mirror  and  the  KDP 
crystal.  From  tabulated  dispersion  data  for 
air,  a  phase  shift  of  180°  occurs  for  an  air 
pressure  equivalent  to  22  cm  Hg.  The  second- 
harmonic  intensity  emerging  from  the  box  is 
shown  in  Fig.  2  as  a  function  of  air  pressure. 

It  has  the  expected  periodicity. 

The  KDP  platelet  is  a  (100)  cut,  and  can  be 
tilted  around  the  (001)  direction.  The  polar¬ 
ization  of  the  laser  beam  is  almost  parallel 
to  the  (010)  direction.  In  this  manner  the  sec¬ 


ond-harmonic  production  in  the  platelet  may 
be  kept  very  small  and  made  comparable  to, 
but  slightly  larger  than,  the  second-harmonic 
production  in  the  mirror. 

The  nonlinear  mirror  may  be  turned  around 
its  (110)  normal  so  that  the  laser  field  points 
along  the  (001)  cubic  axis.  In  this  case  no  sec¬ 
ond  harmonic  is  generated  in  the  mirror.4 
The  KDP  platelet  is  then  tilted,  so  that  its 
second-harmonic  output  is  a  relative  maximum. 

Next  the  mi’Tor  is  turned  so  that  the  funda¬ 
mental  field  points  along  the  (110)  direction. 

The  mirror  then  produces  a  second- harmonic 
field  parallel  to  that  generated  in  the  KDP  crys¬ 
tal.  The  interference  pattern  in  Fig.  2  is  ob¬ 
served,  when  the  air  pressure  in  the  cell  is 
increased.  If  the  fundamental  field  points  along 
the  (111)  direction  o.r  the  crystal  mirror,  no 


FIG.  2.  The  Interference  between  the  second-harmonic  fields  produced  by  the  same  laser  beam  in  the  GaAs 
mirror  and  the  KDP  platelet,  as  a  function  of  the  air  pressure  in  the  box  shown  in  Fig.  1. 
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interference  is  observed  because  in  this  case 
the  fields  generated  in  the  two  samples  are 
orthogonal. 

The  position  of  the  first  extremum  in  Fig.  2 
is  a  measure  for  the  phase  relationship  between 
the  squared  fundamental  field  and  the  second- 
harmonic  field  after  reflection  from  the  non¬ 
linear  mirror.  From  the  angles  for  which  ze¬ 
ro  harmonic  production  occurs  in  KDP,  one 
can  determine  exactly  the  phase  shift  ftV’KDP 
of  the  harmonic  generated  in  KDP  with  respect 
to  that  of  the  harmonic  from  the  nonlinear  mir¬ 
ror.  For  our  particular  case  6^KDP=^  +13°30') 
±  3°30'  when  the  KDP  is  tilted  so  that  its  sec¬ 
ond-harmonic  output  is  a  relative  maximum. 

The  sign  of  the  interference  is  reversed  if 
either  the  33 m  crystal  or  the  KDP  crystal  is 
replaced  by  its  piezoelectric  antipode,  or  in¬ 
version  image.  All  phases  in  the  following  dis¬ 
cussion  have  therefore  an  ambiguity  of  180°. 

This  mnbiguity  of  sign  can  only  be  eliminated 
if  the  sign  of  the  piezoelectric  tensor  is  mea¬ 
sured  and  related  to  the  absolute  structure, 
both  for  KDP  and  the  43m  crystals. 

The  experimental  phase  shift  A  between  the 
position  for  zero  pressure  and  the  first  ex¬ 
tremum  can  be  determined  from  Fig.  2  to  be 
-123°±  4°(+180°)  for  GaAs.  For  a  mirror  of 
InAs  this  phase  angle  was  found  to  be  -90° 

±  4°(+180°)  and  for  ZnTe  0°*  4°(+180°).  From 
these  observations  the  phase  angle  of  the 
complex  nonlinear  susceptibility  may  be  cal¬ 
culated  from  the  relation 

A=V,NL  +  V,(2u;)“2^(a!)“^KDP'  ^ 


The  phase  angle  <p(u>)  is  the  phase  shift  on  re¬ 
flection  of  the  fundamental  field,  determined 
by  the  Fresnel  formula  for  45°  angle  of  inci¬ 
dence, 


E  (w)  1— Qf  (or) 

_ 

EX  w)  1  +a(<A>) 


iifi'jj) 


(2) 


where 


a(o))=;{2e(w)-l}1/2. 


The  phase  angle  ^(2a<)  follows  from  the  non¬ 
linear  laws  of  reflection,®  which  for  the  geom¬ 
etry  of  Fig.  T  may  be  written  in  the  form 


Er (2w)  =  E.2( uit  I  \NL  1  expOVNL)>-2£’^  ^ 2lJ}\ 


with 


r^e 


i<p(  2w) 


_  3277  _ 0(20.0 _ _ _  .  . 

"v2~  {a(uO  +  l}2{a(2o.O  +  l}2{o(u.O  +a(2w)}  ' 

When  the  known  values  of  the  complex  dielec¬ 
tric  constants  are  inserted  in  Eqs.  (l)-(3),  the 
following  phase  angles  for  the  nonlinear  sus¬ 
ceptibility  are  found:  <£nl  =  -71°±  19°(+180°) 
for  GaAs,  <?NL  =  -37°±  20°(+180°)  for  InAs, 
and  <^NL=  34°i  24°(+180°)  for  ZnTe.  With  the 
previously  measured  second-harmonic  inten¬ 
sities2'3  the  estimated  nonlinear  susceptibili¬ 
ties  relative  to  KDP  become 


NL 


lGaAs 


/xKDP  =  ±(165“475i)’ 


NT 

Ws  -W*(272-2Mi)’ 

‘znTe^KDP'1''552*372'’' 


♦This  research  was  supported  by  the  Joint  Services 
Electronics  Program  under  contract  Nonr  I860  (16). 

tOn  leave  from  Harvard  University. 

*N.  Bloembergen,  Proe.  IEEE  151,  124  (1963) 

2N.  Bloembergen,  R.  K.  Chang,  J.  Ducuing,  and 
P.  Lallemand,  Proceedings  of  the  International  Con¬ 
ference  on  Semiconductors,  Paris,  1964  (Academic 
Press,  Inc.,  New  York,  1964). 

3R.  K.  Chang,  J.  Ducuing,  and  N.  Bloembergen  (un¬ 
published). 

4J.  Ducuing  and  N.  Bloembergen,  Phye.  Rev.  Letters 
10,  474  (1963). 

®N.  Bloembergen  and  P.  S.  Pershan,  Phys.  Rev.  128, 
606  (1962).  Especially  Eq.  4.13  of  this  paper  is  used 
with  a  +  =  jt/2  and  6 ^  =  n/4. 
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DISPERSION  OF  THE  OPTICAL  NONLINEARITY  IN  SEMICONDUCTORS* 

R.  K.  Chang 

Gordon  McKay  Laboratory.  Harvard  Univcri.ty,  Cambridge,  Massachusetts 

and 

J.  Ducuing 

Physics  Department,  Massa^  husetts  Institute  of  Technology,  Cambridge,  Massachusetts 

and 

N.  Bloemberger.t 

Department  of  Electrical  Engineering  and  Physics,  University  of  California,  Berkeley,  California 

(Received  8  June  1965) 


The  dispersion  of  the  modulus  of  the  nonlin¬ 
ear  susceptibility  (where  NL  stands 

for  nonlinear)  responsible  for  second-harmon¬ 
ic  (SH)  generation  has  been  directly  measured 
for  the  first  time.  The  crystals  investigated 
were  several  semiconductors  having  zincblende 
(43w)  symmetry.  The  frequency  range  not  on¬ 
ly  includes  the  absorption-edge  region,  but  al¬ 
so  extends  well  inside  this  region.  A  rapid  vari¬ 
ation  of  lxu^k(2u>)l  as  a  function  of  frequen¬ 
cy  is  observed.  Pronounced  maxima  can  be 
correlated  with  the  presence  of  the  critical 
points  in  the  joint  density  of  states  in  the  vi¬ 
cinity  of  the  fundamental  and/or  harmonic  pho¬ 
ton  energy.  Such  correlations  are  similar  to 
the  ones  found  in  the  linear  optical  data.1'2  In 
addition,  the  nonlinear  susceptibility,  which 
is  zero  when  the  crystal  has  inversion  symme¬ 
try,  depends  on  the  antisymmetric  potential3 
which  mixes  the  wave  functions  of  different 
parity.  This  then  leads  to  marked  differences 
between  the  linear  and  nonlinear  cases. 

Nine  different  fundamental  sources  having 
photon  energies  in  the  range  from  1.17  to  2.34 
eV  were  used.  The  lowest  photon  energy  is 
that  of  a  Nd-doped  glass  laser,  the  highest  that 
of  its  harmonic  obtained  by  prior  frequency 
doubling  in  a  KDP  crystal.  The  other  frequen¬ 
cies  are  those  of  a  ruby  laser  and  its  Stokes 
lines  obtained  by  stimulated  scattering  from 
various  fluids:  H2  gas  gives  a  Stokes  energy 
of  1.27  eV,  cyclohexane  1.43  eV:  oenzene  gives 
a  first  Stokes  energy  at  1.66  eV  and  a  second 
Stokes  at  1.54  vV;  nitrobenzene  gives  a  Stokes 
line  at  1.62  eV;  and  CS2  at  1.70  eV. 

The  experimental  arrangement  is  shown  in 
Fig.  1.  The  nonlinear  susceptibility  of  single 
crystals  was  measured  by  the  method  of  re¬ 
flected  harmonics.4  Second-harmonic  produc¬ 
tion  in  a  quartz  plate  is  us^d  to  monitor  the 


□ 


PHOTOMULTIPLIER  EM!  6J55B 

interference  filter 


FIG.  1.  Diagram  of  the  experimental  arrangement 
for  the  case  when  the  incident  radiation  Is  produced 
by  stimulated  Raman  scattering. 


laser  intensity.*  The  intensity  of  the  reflected 
harmonic  beam  from  the  semiconductor  is  mea¬ 
sured  and  compared  to  that  produced  In  trans¬ 
mission  through  a  KDP  crystal  by  the  same 
source.*  The  nonlinear  susceptibility  of  KDP 
is  assumed  to  be  independent  of  frequency  in 
this  range  and  taken  as  a  reference. 

The  nonlinear  second-order  sysceptibility 
tensor  for  a  crystal  with  43w  symmetry  has 
only  three  elements  diffe-ent  from  zero: 
XuNL(2w)  =  x*sNL(2u>)  =  X3#Nk(2u>).  An  expres¬ 
sion  for  the  nonlinearity  appropriate  for  a  semi¬ 
conductor  has  been  derived.7**  A  straightfor¬ 
ward  rearrangement  of  this  expression  yields 


NL 


*14 


(2ui)  =  A  lim  XL 

Imw-0  1,1' 


d3kf(E .  ) 
BZ 


wll'  *  ' 

where  A  is  frequency  independent, 
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interference  is  observed  because  in  this  case 
the  fields  generated  in  the  two  samples  are 
orthogonal. 

The  position  of  the  first  extremum  in  Fig.  2 
is  a  measure  for  the  phase  relationship  between 
the  squared  fundan  ental  field  and  the  second- 
harmonic  field  after  reflection  from  the  non¬ 
linear  mirror.  From  the  angles  for  which  ze¬ 
ro  harmonic  production  occurs  in  KDP,  one 
can  dete.  i  ine  exactly  the  phase  shift  S^KDP 
of  the  harmonic  generated  in  KDP  with  respect 
to  that  of  the  harmonic  from  the  nonlinear  mir¬ 
ror.  For  our  particular  case  5v?kdp  ~  fa  +  13°30') 
±  3°30'  when  the  KDP  is  tilted  so  that  its  sec¬ 
ond-harmonic  output  is  a  relative  maximum.  • 

The  sign  of  the  interference  is  reversed  if 
either  the  43mi  crystal  or  the  KDP  crystal  is 
replaced  by  its  piezoelectric  antipode,  or  in¬ 
version  image.  All  phases  in  the  following  dis¬ 
cussion  have  therefore  an  ambiguity  of  180°. 

This  ambiguity  of  sign  can  only  be  eliminated 
if  the  sign  of  the  piezoelectric  tensor  is  mea¬ 
sured  and  related  to  the  absolute  structure, 
both  for  KDP  and  the  ?3»;i  crystals. 

The  experimental  phase  shift  A  between  the 
position  for  zero  pressure  and  the  first  ex¬ 
tremum  can  be  determined  from  Fig,  2  to  be 
-123°±  4°(+180°)  for  GaAs.  For  a  mirror  of 
InAs  this  phase  angle  was  found  to  be  -90° 
t  4°(+l?.0°)  and  for  ZnTe  0°±  4°(+180°).  From 
these  observations  the  phase  angle  <?NL  of  the 
complex  nonlinear  susceptibility  may  be  cal¬ 
culated  from  the  relation 


The  phase  angle  is  the  phase  shift  on  re¬ 
flection  of  the  fundamental  field,  determined 
by  the  Fresnel  formula  for  15°  angle  of  inci¬ 
dence, 


Er (w)  l-a(w) 
E.(w)  l+a(u>) 


■r,e 


(2) 


where 


a(w)  =  {2e(u.<)-l},/2. 


The  phase  angle  <p( 2u>)  follows  from  the  non¬ 
linear  laws  of  reflection,*  which  for  the  geom 
etry  of  Fig.  1  may  be  written  in  the  form 


£/J(2u.)  =  Ej.2(w)lxNL!exp(/VNL)r2c,V  (2u,), 


with 


i<p(  2u>) 


32rr _ Of(2a>) _ 

=  7T  {o(ce)  +  l[2{a(2*)  +  l}2{a(u>)  +  a(2a>)}‘ 

When  the  known  values  of  the  complex  dielec¬ 
tric  constants  are  inserted  in  Eqs.  (l)-(3),  the 
following  phase  angles  for  the  nonlinear  sus¬ 
ceptibility  are  found:  </j^L  =  -7ra  19°(+180°) 
for  GaAs,  ^NL  =  -37°±  20°(+180°)  for  InAs, 
and  </>nL  =  34°±  24°(+180°)  for  ZnTe.  With  the 
previously  measured  second- harmonic  inter.- 
sities2,i  the  estimated  nonlinear  susceptibili¬ 
ties  relative  to  KDP  become 


NT 

X„  .  /x_  =  ±(165-475z), 

AGaAs  AKDP  * 


NI 

‘.nAs  AKdp=‘«’2-2°«). 
W*KDP=*<552  *372i)’ 


♦This  research  was  supported  by  the  Joint  Services 
Electronics  Program  under  contract  Nonr  1SG6  (1G). 

tOn  leave  from  Harvard  University. 
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The  dispersion  of  the  modulus  of  the  nonlin¬ 
ear  susceptibility  xu^'(2w)  (where  NL  stands 
for  nonlinear)  responsible  for  second- harmon¬ 
ic  (SH)  generation  has  been  directly  measured 
for  the  first  time.  The  crystals  investigated 
were  several  semiconductors  having  zincblende 
@3»i)  symmetry.  The  frequency  range  not  on¬ 
ly  includes  the  absorption-edge  region,  but  al¬ 
so  extends  well  inside  this  region.  A  rapid  vari¬ 
ation  of  lxu^k(2w)l  as  a  function  of  frequen¬ 
cy  is  observed.  Pronounced  maxima  can  b  ; 
correlated  with  the  presence  of  the  critical 
points  in  the  joint  density  of  states  in  the  vi¬ 
cinity  of  the  fundamental  and/or  harmonic  pho¬ 
ton  energy.  Such  correlations  are  similar  to 
the  ones  found  in  the  linear  optical  data.1'2  In 
addition,  the  nonlinear  susceptibility,  which 
is  zero  when  the  crystal  has  inversion  symme¬ 
try,  depends  on  the  antisymmetric  potential1 
which  mixes  the  wave  functions  of  different 
parity.  This  then  leads  to  marked  differences 
between  the  linear  and  nonlinear  cases. 

Nine  different  fundamental  sources  having 
photon  energies  in  the  range  from  1.17  t  2.34 
eV  were  used.  The  lowest  photon  energy  is 
that  of  a  Nd-doped  glass  laser,  the  highest  that 
of  its  harmonic  obtained  by  prior  frequency 
doubling  in  a  KDP  crystal.  The  other  frequen¬ 
cies  are  those  of  a  ruby  laser  and  its  Stokes 
lines  obtained  by  stimulated  scattering  from 
various  fluids:  Ha  gas  gives  a  Stokes  energy 
of  1.27  eV,  cyclohexane  1.43  eV;  benzene  gives 
a  first  Stokes  energy  at  1.66  eV  and  a  second 
Stokes  at  1.54  vV;  nitrobenzene  gives  a  Stokes 
line  at  1.62  eV;  and  CS2  at  1.70  eV. 

The  experimental  arrangement  is  shown  in 
Fig.  1.  The  nonlinear  susceptibility  of  single 
crystals  was  measured  by  the  method  of  re¬ 
flected  harmonics.4  Second -harmonic  produc¬ 
tion  in  a  quartz  plate  is  used  to  monitor  the 


RmO'OWUITipuER  £Ml  6HJ8 
interference  filter 


FIG.  1.  Diagram  of  the  experimental  arrangement 
for  the  case  when  the  incident  radiation  is  produced 
by  stimulated  Raman  scattering. 


laser  intensity.*  The  intensity  of  the  reflected 
harmonic  beam  from  the  semiconductor  is  mea¬ 
sured  and  compared  to  that  produced  in  trans¬ 
mission  through  a  KDP  crystal  by  the  same 
source.*  The  nonlinear  susceptibility  of  KDP 
is  assumed  to  be  independent  of  frequency  in 
this  range  and  taken  as  a  reference. 

The  nonlinear  second-order  sysceptibility 
tensor  for  a  crystal  with  43m  symmetry  has 
only  three  elements  different  from  zero: 
XMNL(2w)  =  X2j^k(2w)  =  Xjs^k(2u>).  An  expres¬ 
sion  for  the  nonlinearity  appropriate  for  a  semi¬ 
conductor  has  been  derived.7'*  A  straightfor¬ 
ward  rearrangement  of  this  expression  yields 

*,  .NL(2tu)=A  ltm  2  J  A/<£  ,) 


(«„.“<*)  V*(wl 

xl"’-“„,’+4"’-V''  (1) 


where  A  is  frequency  independent,  /(Ej^)  is 
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the  Fermi-Dirac  distribution  function  giving 
the  occupation  number  of  the  state  with  the  en¬ 
ergy  Ei  ff.  The  summation  is  over  band  index¬ 
es,  and  the  integration  extends  over  the  Bril- 
louin  zone  (BZ).  The  quantities  <fy,(/fe)  which 
Involve  a  summation  over  intermediate  states 
are  real  and  independent  of  frequency.  They 
vanish  for  a  material  possessing  inversion  sym¬ 
metry.  The  above  expression  is  to  be  compared 
with  the  corresponding  one  for  th°  linear  sus¬ 
ceptibility 
L,  . 

X  ,(2u))=B  lim 
Imw-0 


eL 
(,!•  BZ 


A/<V 


X— i! - -r-= - (2) 


w 


W 


4u>*- 


v 


2> 


ll' 


where  B  is  also  independent  of  frequency,  the 
Pll'lk)  are  suitable  reduced-momentum  matrix 
elements,*  and  L  stands  for  linear.  The  |P^,(ib)l* 
are  known  to  be  slowly  varying  functions  of  k .* 
Were  this  to  be  also  true  of  the  quantities  Q^,(k), 
the  dispersion  of  x**k(2w)  would  be  almost  en¬ 
tirely  governed  by  the  joint  density  of  states 
between  bands  1,1',  as  is  the  case  for  xL-  If 
the  contribution  from  the  highest  valence  band 
and  the  lowest  conduction  band  is  dominant, 
as  in  the  case  of  x^,  the  nonlinear  suscepti¬ 
bility  could  then  be  epresented  to  an  excellent 
approximation  by  a  linear  combination  of  the 
values  x^(w)  and  xL(2w).  In  general,  such  a 
relation  will  not  hold,  as  the  coefficients 
which  involve  triple  products  of  the  reduced- 
momentum  matrix  elements,  depend  on  the  phase 
of  these  latter  quantities.  The  Q’s  can  vary  ap¬ 
preciably  over  the  Brillouin  zone,  changing  in 
sign  as  well  as  in  magnitude.  Because  of  this, 
the  dispersion  of  xNL  will  not  be  determined 
entirely  by  the  joint  density  of  states.  There 
may  also  appear  additional  structure  due  to  the 
variation  of  the  <fyj,(fc)  over  the  Brillouin  zone. 

The  experimental  error  brackets  stem  from 
two  causes.  First,  there  is  a  consistent  error 
of  ±10%  in  determining  the  ratio  of  intensity 
of  the  reflected  SH  from  the  semiconductor  and 
the  intensity  of  the  SH  from  KDP  by  transmis¬ 
sion.  The  second  source  of  error,  which  for 
many  cases  is  the  larger  of  the  two,  results 
from  not  being  ab!e  to  determine  exactly  the 
real  and  imaginary  part  of  the  linear  dielectric 
constant,  especially  when  one  or  both  of  them 
are  sharply  peaked.  The  Fresnel  factors  for 
nonlinear  reflectors  are  a  sensitive  function 


of  them.  The  values  of  the  linear  dielectric 
constant  for  GaAs,  InAs,  and  InSb  were  obtained 
from  Philipp  and  Ehrenreich.*  For  ZnTe,  a 
Kramers-Kronig  transform  was  performed 
using  the  available  data  on  reflectivity  and  suit¬ 
able  extrapolation  for  the  high-energy  region.10-1* 

The  case  of  ZnTe  [Fig.  2(a)]  provides  a  strik¬ 
ing  illustration  of  the  variation  of  lxn^L(2w)l 
when  the  harmonic  photon  energy  penetrates 
progressively  inside  the  absorption  region. 

Over  the  investigated  range  the  absorption  of 
the  fundamental  is  negligible,  but  that  of  the 
harmonic  increases  by  more  than  an  order  of 
magnitude  to  reach  a  maximum  around  2u>  =  3.6 
eV,  which  corresponds  to  the  peak  called  Ey 
in  semiconductor  terminology  in  the  joint  den¬ 
sity  of  states.  Under  these  circumstances  the 
second  term  inside  the  brackets  in  Eq.  (1)  ought 
to  predominate  entirely;  that  is,  the  density 
of  states  at  2fiw  governs  the  behavior  of 
IXmNL(MI.  Experimentally,  as  2fiw  pene¬ 
trates  inside  the  absorption  band,  IxmN^(2w)I 
increases  progressively.  Soref  and  Moos13 
had  previously  observed  a  similar  behavior 
at  a  fixed  light  frequency  when  the  band  gap 
was  varied  in  (ZnCd)  S  alloys.  However,  as 
they  pointed  out,  their  method  assumes  that 
the  properties  cf  the  seir.iconductor-alloy  se¬ 
ries  are  such  that  varying  the  band  gap  by  al¬ 
loying  is  equivalent  to  shifting  the  applied  la¬ 
ser  frequency. 

A  different  situation  exists  in  InSb  [Fig.  2(b)]. 
Here  both  the  harmonic  and  fundamental  pho¬ 
ton  energy  are  absorbed.  As  the  source  fre¬ 
quency  increases,  the  fundamental  photon  en¬ 
ergy  flu)  gets  closer  to  the  E x  peak,  which  for 
InSb  appears  at  a  much  lower  value  (1.83  eV) 
than  for  ZnTe.  On  the  other  hand,  there  ex¬ 
ists  no  critical  point  of  the  joint  density  of  states 
within  the  interval  covered  by  the  harmonic  pho¬ 
ton  energy,  as  evidenced  by  the  fact  that  x^(2w) 
has  no  pronounced  structure.  This  should  re¬ 
sult  in  a  strong  influence  of  the  first  term  in 
Eq.  (1)  on  the  dispersion  of  u(2w).  A  pro¬ 
nounced  variation  of  IXm^^(2w)I  is  observed 
experimentally,  which  seems  to  be  somewhat 
related  to  that  of  the  real  part  of  the  linear  sus¬ 
ceptibility  x^w).  The  maximum  value  appears 
at  fiw  =  1.6  eV,  where  the  Ixh^L(2u>)!  =  5x10-# 
esu.  This  is  the  highest  value  measured  for  a 
second-order  susceptibility,  nearly  2000  times 
higher  than  the  corresponding  value  for  KDP. 

For  GaAs  [Fig.  2(c)],  the  E,  peak  again  ap¬ 
pears  at  higher  value  (2.90  eV),  and  the  main 
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FIG.  2.  The  modulus  IxhNL<2w)|  of  the  nonlinear  susceptibility  of  (a)  ZnTe;  (b)  InSb;  (c)  GaAs;  and  (d)  InAs. 
The  bar  indicates  the  maximum  uncertainty  resulting  from  various  sources  of  error.  For  the  sake  of  comparison, 
the  real  and  imaginary  parts  of  the  linear  susceptibility  xL(2w)  have  been  plotted  in  all  cases,  while  those  of  *L(u>) 
are  also  displayed  in  case  (b). 


features  of  the  dispersion  are  given  by  the  sec¬ 
ond  term  in  Eq.  (1).  Here  a  pronounced  peak 
in  |Xmnl(2w)I  also  exists  in  the  neighborhood 
of  the  Ej  peak.  Two  novel  features  appear  in 
the  results  for  GaAs.  First,  at  the  highest 
source  frequency,  the  harmonic  photon  energy 
is  near  the  Ea  peak  \.AT 5^  —  X 3C»-^ ic)-  A  corre¬ 
sponding  increase  in  the  nonlinearity  is  observed. 
The  contribution  from  the  first  term  of  Eq.  (1) 
might  also  be  significant  at  this  point.  Second, 
a  very  sharp  decrease  in  lxu^^(2w)l  is  observed 
around  1.67  eV.  As  seen  from  the  linear  data, 
the  joint  density  of  states  does  not  exhibit  such 
a  fast  variation  and  thus  cannot  be  held  respon¬ 
sible  for  this  sharp  decrease.  Such  data  must 


probably  be  accounted  for  by  a  pronounced  struc¬ 
ture  in  the  Q’ s.  InAs  [Fig.  2(d)]  yields,  as  ex¬ 
pected,  results  very  similar  to  those  obtained 
for  GaAs.  Both  these  features  were  observed 
again. 

Our  dispersion  results  clearly  show  the  in¬ 
fluence  of  the  critical  points  in  the  density  of 
states.  They  also  bring  out  some  features  re¬ 
lated  to  the  Q’s  and  hence  to  the  hybridization 
of  the  wave  function.  The  dispersion  of  the  real 
and  imaginary14  part  of  lxNL(2u>)l  should  pro¬ 
vide  some  information  not  available  from  the 
study  of  linear  properties  alone. 
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Although  a  great  (leal  of  experimental  work  has 
been  done  in  the  field  of  stimulated  Raman  emission, 
quantitative  measurements  on  the  Raman  gain  have 
been  sparse  The  interpretation  of  such  data3-5  is 
complicated  by  the  multimode  structure  of  the  laser 
beam.6  1  n  this  1  .etter  experimental  evidence  for  these 
multimode  effects  is  presented.  It  is  found  that  the 
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gain  per  unit  length  in  a  Raman  cell  may  vary  by  a 
factor  five  or  more.  The  behavior  of  Raman  oscil¬ 
lators  is  also  profoundly  affected  by  the  mode  struc¬ 
ture  of  the  exciting  laser  beam.  The  unfocused  beam 
of  2-nnn  diani  used  in  our  experiments  has  an 
"ctend'ic",  or  number  of  spatial  modes,  equal  to  25. 
This  number  is  defined  as  the  product  of  the  area  of 
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the  beam  with  the  solid  angle  it  subtends,  divided  by 
the  square  of  the  wavelength.  It  is  important  to  at¬ 
tenuate  the  laser  beam  without  changing  the  mode 
structure.  Since  the  laser  beam  is  linearly  polarized, 
this  can  be  accomplished  by  the  use  of  a  Glan- 
Thompson  prism.  The  beam  passes  through  a 
Raman  oscillator  cell  of  variable  length,  as  one 
window  moved  as  a  piston  in  the  cylindrical  cell.  The 
windows  were  not  coated  and  were  not  optically 
parallel.  No  provision  was  made  to  control  the  feed¬ 
back.  The  emitted  Stokes  intensity  was  measured  by 
a  photocell  with  suitable  filters,  which  accepted  the 
relatively  large  solid  angle  of  6  X  10-3  steradian.  The 
stimulated  Stokes  radiation  is  emitted  in  a  solid 
angle  ofabout2x  10“5 steradian.  Second- and  third- 
order  Stokes  intensities  could  also  be  measured  with 
appropriate  filters.  Figure  1  shows  the  intensity  of 
first  and  second  Stokes  light  in  nitrobenzene  as  a 
function  of  the  length  of  the  Raman  cell  at  constant 
laser  intensity  <5  ti.  The  curves  as  a  function  of  laser 
intensity  at  constant  length  are  essentially  the  same. 


Fig.  1.  Intensity  of  the  first  and  second  Slokcs  lines  of  nitro¬ 
benzene  as  a  function  of  cell  length,  at  constant  laser  intensity. 
The  residual  signal  for  /  =  0  is  due  to  the  finile  transmission  of 
the  fillers  at  the  ruby  frequency.  The  curves  of  the  Stokes  in¬ 
tensity  as  a  function  of  laser  intensity  at  constant  length  l  are 
essential ly  the  same. 


Similar  curves  have  also  been  obtained  for  benzene 
and  other  fluids. 

At  very  low  intensity  only  spontaneous  Raman 
emission  is  observed,  which  is  initially  proportional 
to  3  L.  The  sharp  break  indicates  the  onset  of  re¬ 
generative  feedback.  At  this  point  the  power  gain  in 
a  single  passage  may  be  about  104  to  10s  and  diffuse 
reflection  from  the  windows  is  sufficient  to  obtain 
tegenerative  oscillator  action  in  one  spatial  mode.  If 
only  one  mode  is  oscillating,  the  Stokes  intensity 
should  be  given  by  a  linear  relationship,  A A’  =  C( & l~~ 
d  threshold)  -  This  is  not  the  case.  Stokes  radiation 
conies  out  in  many  modes  at  higher  power  levels, 
which  utilize  more  and  more  modes  or  filaments  of 
the  incident  laser  beam.  The  total  laser  power  is 
never  significantly  depleted.  In  the  case  of  nitro¬ 
benzene  a  small  depiction  is  observed  which  is, 
however,  due  to  two-quanta  absorption.  Apparently 
a  large  fraction  of  the  laser  flux  does  not  participate 
in  Stokes  pumping.  The  observed  Stokes  intensity 
is  well  represented  by 

&S  ~  2  Ct(dL  —  A  (,  threshold) 

There  is  a  di1  iribution  of  thresholds.  Many  modes 
have  a  threohold  three  or  more  times  higher  than 
the  lowest.  This  is  in  agreement  with  the  consider¬ 
ations  of  Bloembergen  and  Shen.6 

The  second  Stokes  intensity  builds  up  in  a  similar 
manner  as  the  Stokes  intensity.  This  second  Stokes 
is  not  generated  parametrically  by  a  polarization 
Pss(vss)  ~  Xst  In  this  case  the  second  Stokes 

radiation  should  come  out  preferentially  in  a  cone 
ittd  ih  UUemUy  t  S3  idwiuk!  be  w  -  / 

This  is  not  observed.  The  second  Stokes  ring  has  an 
intensity  of  the  order  of  1%  of  the  forward  emission 
and  can  be  photographed  only  when  using  a  long 
cell.  7'he  second  Stokes  is  created  by  a  Raman  proc¬ 
ess,  in  which  <3  s  serves  as  pump.  As  soon  as  the 
threshold  is  exceeded  in  the  first  filament,  its  Stokes 
intensity  builds  up  rapidly  depleting  the  laser  power 
in  this  filament.7  It  is  soon  capable  of  creating  3 
in  the  same  filament,  while  other  Stokes  modes  just 
above  their  threshold  are  not  yet  capable  of  doing 
so.  As  the  gain  in  the  Raman  cell  is  increased,  ever 
more  Stokes  modes  arc  fired,  which  in  turn  excite 
an  ever  growing  number  of  second  Stokes  modes, 
etc.  This  is  the  reason  why  the  higher  order  Stokes 
intensities  all  appear  to  build  up  in  a  similar  manner. 
This  effect  has  also  been  observed  by  Bret,'1  and 
Takutna  and  Jennings.3 
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MULTIMODE  EFFECTS  IN  THE  GAIN  OF  RAMAN  AMPLIFIERS 
AND  OSCILLATORS'  II.  AMPLIFIERS 


(Stokes  pumping;  ruby  laser  excitation;  E/T) 


The  high  exponential  gain  and  feedback  in 
Raman  oscillators3  make  a  quantitative  interpreta¬ 
tion  of  experimental  results  obtained  by  their  use 
quite  difficult.  When  the  Raman  gain  is  measured 
instead  with  amplifier  cells,  this  difficulty  disappears 
and  consequently  more  interesting  conclusions  can 
be  drawn  than  was  the  case  for  oscillators.  The 
laser  and  Stokes  radiation  created  in  a  Raman  os¬ 
cillator  traverse  a  variable  distance  d  before  enter¬ 
ing  an  amplifier  cell  of  variable  length  /.  The  laser 
intensity  is  adjusted  so  that  no  regenerative  oscilla¬ 
tion  occurs  in  the  second  cell.  The  Stokes  intensity  is 
attenuated  by  suitable  filters  and  is  measured  at  the 
input  by  a  beam  splitter  and  at  the  output.  The 
aperture  is  confined  so  that  only  Stokes  radiation 
within  the  confines  of  the  laser  beam  is  accepted. 
The  gain  in  the  amplifier  has  been  measured  as  a 
function  of  the  length  /  as  well  as  of  the  separation 
d.  Figure  1  shows  that  the  overall  gain  is  a  sensitive 
function  of  the  separation  d.  Furthermore,  the  gain 
per  unit  length  {d[ln(3.sl£0)]}ldx  in  the  cell  is  con¬ 
siderably  smaller  in  the  first  few  centimeters  than  at 
larger  distances  from  the  cell  from. 

Both  remarkable  facts  are  readily  explained  by 
the  multimode  theory.  If  the  separation  d  between 
the  cells  is  sufficiently  long,  the  intensity  pattern  of 
the  Stokes  radiation  at  the  input  of  the  second  cell 
will  not  have  an  optimum  Correspondence  with  the 
intensity  pattern  at  the  laser  frequency.  A  decom¬ 
position  of  the  Stokes  input  into  new  normal  modes 
corresponding  to  the  laser  intensity  distribution  in 
the  amplifier  cel!  is  necessary.  Many  modes  will  grow 
with  the  average  gain,  which  is  represented  by  the 
initial  slope  of  the  gain  curve.  F.ventually  the  mode 
with  highest  gain  will  take  over.  This  gain  is  rep- 
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resented  by  the  steep  part  of  the  curves.  Since  the 
insertion  loss  for  the  mode  with  maximum  gain  is 


Fig.  1.  Stokes  gain  in  a  second  Raman  cell  in  scries  with  a 
Raman  oscillator  cell  as  a  function  of  amplifier  cell  length  / 
for  different  values  of  the  spacing  d  between  the  two  cells  and 
of  the  input  Stokes  intensity  v 


212 


t 


Volume  6,  Number  10 


APPLIED  PHYSICS  LETTERS 


May  15,  1965 


larger  for  larger  separations  d ,  it  takes  longer  for 
this  mode  to  emerge  for  larger/.  Ordinarily  the  gain 
per  unit  length  is  taken  to  be  proportional  to  the 
Raman  susceptibility  and  the  average  laser  intensity. 
This  gain  corresponds  more  nearly  to  the  initial 
slope  rather  than  the  maximum  slope  in  the  ampli¬ 
fier  cell.  Large  systematic  errors  may  therefore 
have  occurred  in  previous  determinations  of  the 
Raman  susceptibility.  The  effect  of  saturation  due 
to  the  depletion  of  the  laser  power  in  the  mode 
with  maximum  gain  is  evident  in  the  curve  of  Fig.  1 
for  the  larger  Stokes  input  power.  This  saturation 
mechanism  is,  of  course,  responsible  for  limiting 
the  Stokes  power  in  a  Raman  oscillator.  These 
effects  may  be  caused  both  by  the  spatial  and  the 
temporal  intensity  structure  of  the  laser  beam. 

The  spatial  multimode  effect  can  be  demon- . 
strated4  in  a  striking  manner  by  splitting  a  laser 
beam  into  beams  of  equal  intensity  which  excite 
Stokes  radiation  in  identical  Raman  cells,  as  shown 
in  Fig.  2.  If  the  beam  splitter  is  of  good  optical 
quality,  with  a  tolerance  less  than  X/ 20,  the  inten¬ 
sities  of  the  two  Raman  cells  in  parallel  have  a  very 
good  one-to-one  correspondence  in  each  pulse.  The 
relative  rms  deviation  in  the  intensity  ratio &aI&b  is 
about  0.2.  When  the  beam  splitter  has  a  tolerance 
in  flatness  of  X/2,  this  rms  deviation  increases  by  a 
factor  three  to  0.6.  In  this  case  significant  differ¬ 
ences  in  the  mode  structure  of  the  split  beam  in  the 
two  cells  are  present  during  each  pulse. 

The  Raman  amplifier  cell  will  be  very  useful  fora 
quantitative  study  of  the  gain  as  a  function  of  the 
state  of  polarization  of  the  Stokes  and  laser  polar¬ 
ization.  In  principle  all  three  constants5  of  the 
Raman  susceptibility  tensor  iri  isotropic  fluids  may 
be  measured.  It  may  also  be  possible  to  measure  the 
Raman  gain  at  Stokes  frequencies  which  are  not 
emitted  in  Raman  lasers,  because  the  corresponding 


Fig.  2.  Experimental  set-up  to  show  the  influence  of  the 
spatial  multimode  structure  of  the  laser  light  on  Stokes  produc¬ 
tion.  The  diagrammatic  insert  shows  the  scattering  in  the  ratio 
of  the  Stokes  intensities  in  cells  A  and  B  from  each  laser  pulse. 

vibrations  give  too  high  a  threshold.  The  amplifier 
cell  will  also  be  used  by  us  to  unravel  further  the 
relationships  in  mixtures  of  Raman  fluids,  which 
have  recently  been  reported  by  Kaiser,  Maicr,  and 
Giordmaine.6  When  a  more  powerful  laser  beam 
with  a  larger  cross  section  becomes  available,  the 
angular  dependence  of  the  gain  at  Stokes  and  anti- 
Stokes  frequencies7  may  also  be  studied  more 
quantitatively  with  the  amplifier  cell. 
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Light  Waves  with  Exponential  Gain® 
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The  coupled-wave  formalism  is  reviewed  to  give  a  unified  de¬ 
scription  of  parametric  down  conversion  of  light,  the  stimulated 
Brillouin  and  Raman  effects,  stimulated  light  scattering  in  plasmas, 
etc.  The  exponential  character  of  the  gain  is  emphasized.  The  feed¬ 
back  in  oscillators  and  the  filamentary  or  multimode  structure  of  the 
light  beams  present  serious  difficulties  for  quantitative  measure¬ 
ments,  Some  experimental  data  on  the  gain  in  a  Stokes  amplifier 
cell  are  presented.  The  exponential  gain  constant  for  the  benzene 
Stokes  line  has  been  measured  as  a  function  of  polarization  and  of 
concentration  in  liquid  mixtures.  The  relationship  to  the  spontaneous- 
emission  cross  section  in  these  eases  is  discussed. 


THREE  COUPLED  LIGHT  WAVES 

The  nonlinearity  of  an  electromagnetic  medium  leads  to  a  coupling 
between  waves,  which  are  the  eigensolutions  of  the  wave  equation  in 
the  absence  of  the  nonlinearity.1  The  coupling  between  three  electro¬ 
magnetic  waves  was  considered  in  the  theory  of  microwave  para¬ 
metric  amplifiers,  and  the  coupling  between  three  and  more  light 

*This  research  was  supported  by  Project  Defender  of  ARPA  and 
the  Office  of  Naval  Research. 
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waves  was  treated  early  in  the  development  of  nonlinear  optics.1-3 
The  three  nonlinear  coupled  equations  in  a  piezoelectric  crystal  can 
be  written  in  the  form 


v*e\  .  E*  -  -  x*nl  E',  Eg  ciAk-  r 

c  c 

V*EL+  -^T^EL  -  XNLEiEs 


(1) 


Often  the  laser  amplitude  El  may  be  considered  as  a  constant 
parameter,  since  the  reaction  on  it  by  the  much  weaker  fields  at  Es 
and  Ef  is  negligible.  For  Ak  =  0  and  idler  damping  large  compared 
to  the  signal  damping,  one  obtains  a  power  gain  at  u)S 


gs  =  -2<*j 


(4jru>ju.s  \ 


IxnlI?|el|2 

kI>s“i 


(2) 


If  the  exponential  gain  is  sufficient  to  overcome  coupling  and  dif¬ 
fraction  losses,  the  system  will  become  unstable  and  oscillations 
will  start. 

Giordmaine  and  Miller4  have  very  recently  succeeded  in  oper¬ 
ating  a  tunable  light  oscillator  based  on  this  parametric  down- 
conversion  instability.  The  phase-matching  condition  in  a  LiNbG, 
crystal  can  be  realized  for  three  collinear  light  beams.  The  mo¬ 
mentum  matching  depends  sensitively  on  the  temperature  and  allows 
thermal  tuning  of  the  parametric  oscillator. 


COUPLING  BETWEEN  TWO  LIGHT  WAVES  AND  ONE  VIBRA¬ 
TIONAL  WAVE 

The  stimulated  Raman  effect  and  the  stimulated  Brillouin  effect 
can  be  described  in  the  same  manner.1’5  One  must  replace  the  idler- 
wave  equation  in  the  coupled  set  by  the  wave  equation  for  acoustical 
phonons  in  the  case  of  the  Brillouin  effect  and  by  the  wave  equations 
for  optical  phonons  in  the  case  of  the  stimulated  Raman  effect.  The 
coupling  constants  are  now  proportional  to  the  photoelastic  constant 
and  the  electronic -vibrational  coupling,  respectively,  it  is  also  pos¬ 
sible  to  replace  the  idler  wave  by  a  plasma  wave  or  spin  wave.  This 
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coupling  of  light  waves  with  plasmons  and  magnons  is  discussed 
elsewhere  in  these  proceedings. 

The  set  of  coupled  equations  for  the  optical  phonon  wave  Qv  and 
the  Stokes  wave  E8  is 


0V2Q*  +  (4  ~  4)  Qv  +  2iu»vrvQ;  =  N(aa/dQ)ELEs  (3) 

V2E8  +  (€8c4/c2)E8  -  -4tt(u,2 /c2)N(aa/9Q)QjEL  (4) 


with  Qv  =  R(2p)^2,  where  R  is  the  internuclear  distance  and 
p  =  Npv  is  the  reduced  mass  density  of  the  vibration  for  N  mole¬ 
cules  per  unit  volume,  each  with  a  reduced  vibrational  mass  pv- 
The  group  velocity  du.</dkv  of  the  optical  phonons  is  very  small.  The 
frequency  of  the  optical  phonon  is  essentially  independent  of  the  di¬ 
rection  of  the  Stokes  wave,  wv  =  u>o  “0(kL  ~  k8)2/u>o  u>o*  For  for¬ 
ward  Stokes  scattering  the  phase  velocity  of  the  optical  phonon  is 
equal  to  the  phase  velocity  of  the  light  waves.  For  backward  scat¬ 
tering  the  phase  velocity  is  about  twenty  times  ‘•'ower  for  typical 
vibrational  frequencies.  The  steady-state  gain  is  now  always 
reached  during  the  giant  laser  pulse.  Because  of  the  low  group 
velocity  and  the  relatively  strong  vibrational  damping,  the  exponen¬ 
tial  Stoke*  power  gain  per  unit  length  may  be  written 


8s  = 


As 


N2(3a/aQ)2 

2(*iy  F  y 


(5) 


Here  \8  is  the  Stokes  wavelength  in  the  medium.  The  Stokes  Raman 
susceptibility  has  been  introduced  by 


N2(8a/9Q)2  _  N<9o/9R)2 

2a)VFv  2pvu.'vrv 


COUPLING  WITH  ANTI-STOKES  WAVES 

The  vibrational  wave  and  the  laser  wave  also  couple  to  the  anti- 
Stokes  wave  at  u>a  =  *-  uv.  If  initially  only  the  laser  amplitude 

is  different  from  zero,  the  anti-Stokes  waves  cannot  be  generated. 
This  case  of  parametric  up-conversion  has  also  been  discussed  ex¬ 
tensively  for  three  light  waves.  When  the  vibrational  mode  is 
heavily  damped  and  the  solution  analogous  to  the  Stokes  wave  is 
written  down,  one  finds  that  the  sign  of  g8,  or  equivalently  the  sign 
of  the  imaginary  part  of  the  nonlinear  susceptibility,  has  changed. 
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There  is  exponential  loss  at  the  anti -Stokes  frequency,  which  has 
been  culled  stimulated  absorption  by  Stoicheff  and  coworkers.6 

When  the  momentum -matching  conditions  for  Stokes  and  anti- 
Stokes  waves  are  simultaneously  satisfied,  kv  =  kL,  ~  -  ka  ~  kL. 

a  set  of  four  coupled  waves  at  o?l*  wv.  ws»  an<*  ^'a  rnust  be  solved 
simultaneously.  If  the  laser  field  is  again  treated  as  a  constant 
parameter  and  the  vibrational  wave  is  eliminated  because  of  heavy 
damping  and  lo*.  group  velocity,  one  retains  a  set  of  two  linearized 
coupled  equations  for  Es  and  Ea. 


V2ES  +  ts(u>2 /c2)Es  =  4ff(u)2/c2)lXslELl2Es  +  <XaXs>  l/2E2LE*l 
r2Ej(^/c2)Ej  =  47r(o;2/c2)|(xaXs),/2E*L2Es  +  xJIElAe* 


(7) 


There  is  one  mode  with  exponential  gain  which  has  an  appreciable 
ar.ti-Stokes  character  near  the  phase-matched  direction.  This  expo¬ 
nential  gain  at  *ca  explains  the  copious  production  of  anti -Stokes 
rings.  The  formal  solution  gives  zero  exponential  gain  for  exact 
matching  k«?  +  ka  =  2kL  of  the  unperturbed  propagation  vectors  in  the 
medium.7  In  this  case  each  of  the  two  coupled  eigenmodes  has  50  per 
cent  Stokes  and  anti-Stokes  character.  This  point  has  been  discussed 
by  several  authors.*  It  is  a  rather  academic  point,  however,  because 
the  dip  to  zero  has  rn  extremely  narrow  angular  width  AO  «  Ak/k 
~  gsA.  Since  the  gain  gs  rarely  exceeds  unity,  i.e.,  an  e -folding 
length  of  1  cm,  one  has  AO  <  10"4  rad.  A  finite  cross  section  of  the 
beam,  if  it  is  confined  to  less  than  1  cm  in  diameter  as  is  always 
the  case  experimentally,  or  a  nonuniform  intensity  distribution 
across  the  beam  as  would  be  the  case  if  gain  in  lateral  directions 
is  allowed,*  will  obliterate  this  sharp  dip  in  the  gain.  There  is  ex¬ 
ponential  gain  at  the  anti -Stokes  frequency  near  the  phase-matched 
direction,  and  as  a  result  the  Stokes  character  of  the  wave  is  not 
100  per  cent.  There  is  a  drop  in  the  Stokes  intensity  near  the 
phase-matched  direction,  as  has  been  observed  by  Chiao  and 
Stoicheff9  and  by  Garmire.*0 

The  Stok.  s  and  anti -Stokes  waves  can  acquire  such  intensities 
that  they  noticeably  deplete  the  laser  beam  and  can  in  turn  create 
higher-order  waves.  The  set  of  coupled  wave  equations  then  becomes 
very  complex. 


MULTIMODE  EFFECTS 

This  complexity  is  compounded  by  the  fact  that  the  incident  laser 
beam  usually  cannot  be  represented  by  a  single  homogeneous  mono- 
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chromatic  plane  wave.  The  Stokes  and  anti-Stokes  radiation  of  var¬ 
ious  orders  is  observed  to  have  a  distribution  of  directions  and  fre¬ 
quencies  that  implies  the  participation  of  a  very  large  number  of 
modes.11  As  in  all  instability  problems,  the  modes  with  the  maxi¬ 
mum  gain  have  to  be  found. 

In  practice  the  laser  and  Stokes  intensity  patterns  cannot  be 
described  simply  as  the  interference  pattern  of  two  or  more  plane 
waves.  As  will  be  discussed  in  more  detail  in  the  concluding  section, 
after  presentation  of  experimental  data,  one  may  expect  patterns  of 
hot  filaments,  changing  with  time  during  the  pulse.  Because  of  the 
nonlinear  nature  of  the  phenomena,  the  observed  gain  for  a  given 
total  laser  output  can  be  an  order  of  magnitude  larger  han  the  gain 
calculated  for  a  uniform  field. 


POLARIZATION  MODKS.  RELATION  BETWEEN  THE  RAMAN 
TENSOR  AND  SPONTANEOUS  EMISSION 

Even  if  the  direction  and  the  frequency  of  the  laser  and  Stokes 
waves  are  completely  defined,  there  remains  the  twofold  degeneracy 
of  the  polarization.  One  therefore  has  strictly  speaking  always  to 
solve  at  least  the  two  coupled  equations  for  the  two  Stokes  polari¬ 
zations.  Consider,  for  simplicity,  the  case  of  a  Raman  liquid.  In  an 
isotropic  medium  the  Stokes  polarization  is  completely  specified  by 
three  components  of  the  Raman  susceptibility  tensor, 


Raman  £  \  .  Raman 

<  . nEL(E,.-Es> 


p(«s)  ES(E|.-E*L)  *  E,(ES-El) 

Raman 


It  should  be  borne  in  mind  that  the  Raman  susceptibility  is  a  complex 
quantity  which  obeys  the  same  Kramers -Kronig  relation  as  a  linear 
susceptibility.  Exactly  at  resonance  it  is  pure  imaginary.  The  real 
part  is  often  equally  important,  since  the  observed  frequency  distri¬ 
bution  of  higher-order  Stokes  and  anti -Stokes  radiation  shows  that 
one  often  works  in  the  wings  of  the  vibrational  resonance.  All  equa¬ 
tions  retain  their  validity  for  complex  x- 

Consider  the  case  that  both  the  laser  and  the  Stokes  wave  propa¬ 
gate  in  the  z  direction.  The  amplitudes  Exl  and  EyL  are  treated  as 
fixed  parameters.  The  set  of  coupled  equations  for  the  two  Stokes 
polarizations  is 
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^2^XS  +  *S^Wg  A^Ejjg 

-(4ffu)g/e2)((xi  f  X2  f  X3) !  ExlI2  +  X  jlKyLl2lExs 
f  (4jro;|/c2)(x2ExLEyL  +  X3EjLEyL)-ys 

^‘Eys  +  (  g(u>g  / C2  )  Eyg  = 

-<4,4/c2)(X2EjLEyL  -  X3E*LExL)Exs 
+  (4jtu)2/c2)|  Xi|ExL|‘j  Mxi  +  X2  +  Xa)|Ey  i^llEyg 

The  eigenmodes  may  be  solved  for  the  case  of  general  elliptical 
polarization  of  the  incident  laser  beam.  Two  particular  cases  are 
evidently  of  special  interest. 

1.  The  laser  beam  is  linearly  polarized. 

In  this  case  we  may  put,  without  loss  of  generality,  Evl  EJl 
3  0.  The  Stokes  beam  Esx  polarized  parallel  to  the  laser  beam 
has  a  power  gain 

8s,ll  =  -(S^sAsM*','  f  X2  +  Xj ) I  ExL|2  .  (8) 

This  is  the  Stokes  radiation  that  has  been  observed  in  Raman  oscil¬ 
lators.  The  perpendicular  eigenmode  %s  has  a  gain 

8s, 1  *  ”  {8jr2t s A«)X  ?  I  ExjJ 2  (9) 

In  a  Raman  oscillator  only  the  linear  mode  with  the  larger  gain  will 
be  excited.  The  corresponding  real  parts  represent  a  birefringence 
of  the  liquid,  induced  by  the  laser  beam.  Terhune13  has  observed  an 
intensity -dependent  index  of  refraction  on  laser  beam  itself.  The 

effect  should,  however,  be  particularly  large  on  a  Stokes  line  just  off 
the  vibrational  resonance,  where  the  real  part  is  enhanced. 

2.  The  laser  beam  is  circularly  polarized. 

In  this  ease  we  may  take  Ex^  real  and  write  Kxl  ~  iEyL,  =  E0l. 
Solution  of  the  coupled  wave  equations  shows  that  the  Stokes  eigen¬ 
modes  are  now  also  circularly  polarized.  The  Stokes  mode  with  the 
same  sense  of  polarization  as  the  laser  beam  has  a  power  gain 

Ss.same  “  ~(87r2tg/Ag)  Eql,( x  1  +  X2)  (10) 

while  the  mode  with  the  opposite  sense  of  polarization  has  the  gain 
eonstjmt 


8s, opposite  ”  <8*2*sAs)  EoL<Xi  *  X3) 


01) 
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The  output  of  a  Raman  oscillator,  pumped  with  a  circularly  polarized 
laser  beam,  will  itself  be  circularly  polarized.  Similar  considera¬ 
tions  hold  for  the  stimulated  Brillouin  effect  in  liquids. 

The  real  parts  of  these  susceptibilities  now  lead  to  the  phenom¬ 
enon  of  light-induced  optical  activity.  In  the  presence  of  circular 
polarized  laser  light,  the  plane  of  polarization  of  n  linear  polarized 
beam  at  another  frequency  will  turn.  This  effect  should  again  be 
pronounced  for  Stokes  light  just  off  the  vibrational  resonance,  where 
the  real  parts  of  these  susceptibilities  will  have  the  largest  differ¬ 
ence.  In  principle  the  three  components  can  be  determined  by  meas¬ 
uring  the  gain  constants  in  these  geometries. 

It  has  been  shown  in  Eq.  (6)  that  the  Raman  susceptibility  is  pro¬ 
portional  to  the  derivative  of  the  electronic  polarizability  of  a  mole¬ 
cule  with  respect  to  internuclear  distance  da/dQ.  In  the  same  man¬ 
ner,  the  tensorial  or  anisotropic  properties  of  this  susceptibility 
are  related  to  the  anisotropic  properties  of  the  polarization  deriva¬ 
tives  of  a  molecule,  suitably  averaged  over  all  molecular  orienta¬ 
tions  in  the  liquid.*4  This  procedure  was  developed  by  Placzek15  and 
has  been  reviewed  in  books  on  spontaneous  Raman  scattering. 16 

Consider  an  anisotropic  molecule  with  principal  polarizabilities 
a,.  a2.  and  a3.  Define  a  scalar  part  a  and  an  anisotropic  part  b  by 

a  =  ^-(<*1  +  o2  +  <*3) 

"  "1  1/2 
b  "  -£(<*|  -a2)2  *  -|(a2-a3)2  +  -|(a3-al)2J 

Take  the  derivatives  of  these  quantities  with  respect  to  the  normal 
coordinate  Q  of  the  totally  symmetric  vibration  of  the  molecule 

a  =  a0  +  (da/dQ)  Q  =  a0  *  a'Q 
b  ~  b0  +  (db/dQ)Q  =  b0  +  VQ 

Calculate  the  square  of  (do/dQ)  for  the  particular  Raman  polariza¬ 
tion  of  interest  for  a  given  molecular  orientation,  and  average  over 
all  orientations  in  the  isotropic  liquid.  In  this  manner  one  finds  that 
the  Raman  gain  for  the  modes  with  a  linear  polarized  laser  beam  is 
proportional  to 

Xl  *  Xi'  *  Xi  “  -N(45a'2  <  4b'!)/2MvWvrv 
Xi  =  -3Mb/2/2pvwvrv 


(12) 
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For  the  circular  oolarized  modes  the  gain  is  proportional  to 


X?  +  X2  =  ~ N(45a'2  ►  b'2)/3MVu;vrv 

xT +  xa  -6Nb/2/^vwvrv  (13) 

The  three  components  of  the  Raman  susceptibility  tensor  can  there¬ 
fore  be  expressed  in  terms  of  two  independent  quantities  a'  and  b'. 
This  is  true  if  and  ws  are  sufficiently  far  removed  from  elec¬ 
tronic  resonances  so  that  the  antisymmetric  part  of  the  susceptibility 
tensor  can  be  ignored. 15  The  relation  of  a'  and  b'  to  the  symmetry 
properties  of  a  molecule  is  well  known.  For  totally  symmetric  vi¬ 
brations  a'  and  b'  can  assume  arbitrary  values,  except  for  symmetric 
molecules  such  as  CC14  belonging  to  a  cubic  class;  in  this  case  b' 

.=  0.  For  antisymmetric  vibrations  a'  must  vanish. 

Note  that  the  ratios  of  the  stimulated  Raman  gain  factors  are  the 
same  as  the  depolarization  and  inversion  factors  for  spontaneous 
Raman  emission  in  the  corresponding  geometries.  This  remains 
tr*-e  if  interactions  between  molecules  and  local  field  corrections  in 
the  liquid  are  taken  into  account.  Precisely  the  same  correction  fac¬ 
tors  apply  in  both  cases.  The  polarization  tensor  a  must  be  replaced 
by  L(wl)  *<*  *  k(a!g)  where  L(w)  »3  the  Lorentz  local  field  tensor.  For 
a  liquid  without  associative  effects  of  the  neighboring  molecules  the 
Lorentz  correction  factor  on  the  Stokes  intensity  is  the  scalar 
quantity  (e(a>s)  +  2)2[e(u;L)  +  2|2/81 .  This  factor  also  applies  to  the 
spontaneous  Raman  emission.  The  acting  laser  field  is  now  the 
field  inside  a  Lorentz  cavity.  For  the  spontaneous  emission  at  us 
the  mode  wave  function  of  the  vacuum  should  be  replaced  by  the 
mode  wave  function  for  the  dielectric  with  the  cavity.  There  is  con¬ 
sequently  a  similar  correction  factor  at  us-  The  ratio  between 
stimulated  emission  and  spontaneous  emission  remains  nph(u>s)> 
where  this  excitation  number  of  the  field  oscillator  is  proportional 
to  the  classical  field  intensity.  The  relationship  between  the  total 
Raman  scattering  cross  section  for  one  particular  mode  of  Stokes 
polarization,  which  radiates  in  a  cos2  6  pattern  into  the  total  solid 
angle  4tt,  and  the  Raman  susceptibility  for  that  same  mode  of  polar¬ 
ization  is 


rvhwLwsts 


3 

4jt 


(14) 
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THE  OUTPUT  INTENSITY  OF  RAMAN  OSCILLATORS 

The  Stokes  (u>s  =  -  u;v)  and  second  Stokes  (wSs  =  wL  “  2a.'v) 

intensities  of  a  Raman  oscillator  cell  as  a  function  of  laser  power  arc- 
shown  in  Fig.  1.*  This  particular  example  is  for  a  10-cm-long  nitro¬ 
benzene  cell.  The  output  as  a  function  of  cell  length  at  constant  laser 
power  shows  a  nearly  identical  behavior  and  has  been  published  else¬ 
where.  17  This  output  behavior  is  representative  of  instabilities  pro¬ 
duced  by  laser  beams.  Similar  curves  have  been  obtained  by  Bret, 18 
while  Brewer18  has  observed  the  same  trends  for  the  Brillouin- 
shifted  lines. 

Although  the  windows  in  the  particular  cells  with  which  Fig.  1 
was  obtained  were  not  parallel,  the  onset  of  regenerative  feedback 
is  evident  bv  the  shaip  break  in  the  curve.  At  lower  power  levels 
spontaneous  emission  is  observed.  The  feedback  is  perhaps  caused 


Fig.  I.  The  Stokes  and  second  Stokes  output  intensity  of  a  10-cm- 
long  nitrobenzene  Hainan  cell,  us  a  function  of  laser  pump  power. 
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by  random-scattering  centers,  or  it  may  be  caused  by  a  transverse 
instability  of  the  laser  beam  discussed  in  the  concluding  section. 

\  The  following  features  show  that  the  theory  of  pumping  by  a  single 

uniform  laser  mode  is  inadequate. 

1.  Second  (and  higher)  order  Stokes  radiation  appears  at  almost 
the  same  threshold  as  the  fundamental  Stokes  light. 

2.  Above  threshold  the  output  of  a  feedback  oscillator  should  be  a 
linear  function  of  pump  power.  If  the  curve  of  Fig.  1  is  plotted  on  a 
linear  scale,  it  has  a  positive  curvature. 

3.  For  Raman  ceils  longer  than  about  20  cm,  the  threshold  be- 

j  comes  essentially  independent  of  the  reflectivity  and  parallelism  of 

the  windows.  For  very  short  cells,  however,  the  feedback  by  reflec¬ 
tion  is  important. 

4.  The  laser  power  often  does  not  get  significantly  depleted.  The 
Stokes  power  in  long  cells  usually  does  not  exceed  1  per  cent  of  the 
laser  power  in  unfocused  beams. 

5.  The  Stokes  generation  is  not  a  uniquely  defined  process  but 
has  some  stochastic  properties. 

This  last  point  is  proved  experimentally  by  splitting  the  laser 
beam  and  pumping  two  identical  Raman  cells  in  parallel.17  There  is 
no  strict  one-to-one  correspondence  between  the  outputs  of  the  two 
cells.  Since  the  rms  deviations  become  larger  when  the  beam  split¬ 
ter  is  of  poor  optical  quality,  this  suggests  a  filamentary  of  multi  - 
mode  structure  of  the  laser  beam  which  varies  from  pulse  to  pulse. 
Stochastic  variations  in  the  laser-beam  intensity  may  occur  both  in 
time  and  in  space. 

The  Stokes  output  shown  in  Fig.  1  can  be  represented  rather  well 
by  a  linear  combination  of  filaments  with  a  distribution  of  laser 
i  «  powers  rangiug  over  an  order  of  magnitude. 


*s  =  ^threshold) 


This  distribution  may  also  explain  the  behavior  of  the  second  Stokes 
intensity.  In  the  filaments  with  the  highest  laser  intensity,  the  Stokes 
intensity  builds  up  rapidly  to  pass  the  threshold  for  second  Stokes 
generation,  while  in  other  regions  the  nonuniform  Stokes  intensity  is 
not  capable  of  doing  so.  Only  a  fraction  of  the  total  laser  power  is 
active  in  the  stimulated  Raman  process  because  the  exponential  gain 
in  regions  of  lower  laser  intensity  is  not  sufficient.  A  quantitative 
determination  of  Raman  gain  and  susceptibilities  becomes  very  dif¬ 
ficult  under  these  circumstances. 

The  difficulties  are  enhanced  if  one  takes  into  account  the  angular 
and  frequency  distribution  of  the  Stokes  radiation.  The  best  roeas- 
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urements  would  be  made  from  the  threshold  in  very  short  cells  with 
a  well-controlled  feedback  by  reflecting  mirrors.  In  this  case  one 
would  obtain  a  measure  of  the  gain  in  the  “hottest”  filament.  The 
output  of  oscillators  is  of  course  always  determined  by  some  non¬ 
linearity.  The  depletion  of  the  laser  intensity  is  the  most  likely 
limiting  factor  in  Raman  oscillators.  For  quantitative  measurements 
and  an  experimental  check  of  the  exponential  gain,  stable  amplifier 
cells  must  be  preferred. 


EXPONENTIAL  GAIN  IN  RAMAN  AMPLIFIER  CELLS 

It  is,  of  course,  possible  to  prevent  the  instability  from  develop¬ 
ing  by  choosing  the  cell  length  sufficiently  small  and/or  the  laser 
pump  intensity  sufficiently  low.  In  this  case  one  has  a  stable  ampli¬ 
fier  with  zero  output  for  zero  Stokes  input,  except  for  very  weak 
spontaneous-emission  noise.  The  experimental  arrangement  to 
measure  the  gain  of  the  stable  amplifier  is  shown  in  Fig.  2.  Laser 


Fig.  2.  Diagram  of  experimental  arrangement  for  the  gain  measure¬ 
ment  in  a  Raman  amplifier  cell. 

1.  -  ruby  laser 
P  -  Gian  Thomson  polarizer 
O  =  Raman  oscillator  cell 

M,  =  selective  mirror,  with  a  transmission  loss  of  30  db  at  u>s  and 
l  db  at  wl 

F,  Fj,  t'-> ,  and  F3  =  selective  filters  at  u: l  and  uis 
Q  -  quartz  retardation  plate  for  selection  of  polarization  at  u>L  and 

s 

A  -  Raman  amplifier  cell  of  variable  length  1 

PS  -  polarization  separator 

P.M|  -  photomultiplier  monitor  of  Stokes  input 

PMi  and  PM3  -  photomultipliers  monitoring  the  two  polarizations  of 
Stokes  output 
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Fig.  3.  The  Stokes  output  to  input  ratio  or  gain  of  a  benzene  Raman 
amplifier  cell  ns  a  function  of  cell  length  for  three  different  laser 
pump  power  levels.  The  ratio  of  the  laser  intensities  for  curves 
B,  and  C  is  as  1:0.57:0.30,  respectively.  The  intercepts  of  the 
straight  high-gain  parts  with  the  intensity  axis  give  the  insertion 
loss  of  the  Stokes  input  into  the  filament  with  highest  gain. 

and  Stokes  radiation  is  generated  in  a  Raman  oscillator.  A  selec¬ 
tive  reflector  attenuates  the  Stokes  radiation  by  at  least  30  db  and 
the  laser  pulse  by  about  1  db.  The  30-db  padding  at  u>s  isolates  the 
amplifier  from  the  oscillator.  Additional  selective  absorbers  at  u>l 
and  u)s  may  be  inserted.  The  input  Stokes  powrr  is  monitored  by  a 
photomultiplier  with  suitable  selective  filters.  The  output  Stokes 
power  after  the  light  beams  have  traversed  an  amplifier  cell  of 
variable  length  is  measured  by  another  photomultiplier,  which  is 
again  supplied  with  calibrated  filters  and  attenuators  to  ensure  op¬ 
eration  at  a  suitable  intensity  level. 

The  variation  of  the  amplifier  gain,  defined  as  IgUt/Ign.  is  plotted 
in  Fig.  3  as  a  function  of  cell  length  1  for  three  different  laser  power 
levels.  Curve  A  was  obtained  for  full  laser  power.  In  curves  B  and 
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C  the  laser  power  was  attenuated  by  selective  filters  to  5?  and 
36  per  cent  of  its  full  value,  respectively.  It  is  possible  to  exclude 
saturation  effects  experimentally.  The  same  ratio  IsUt/l!.n  is  ob¬ 
tained  if  the  input  intensity  is  changed  by  an  order  of  magnitude. 

The  semilogarithmic  plot  shows  that  there  is  a  region  of  exponen¬ 
tial  gain  with  a  gain  constant  per  unit  length  proportional  to  the  laser 
power,  given  by  the  slope  of  the  experimental  plot.  The  remarkable 
feature  is,  however,  that  the  initial  slope  is  much  smaller.  The  gain 
constant  in  the  first  few  centimeters  of  the  amplifier  is  at  least  an 
order  of  magnitude  smaller  than  farther  on.  Furthermore,  the  re¬ 
gion  of  low  gain  becomes  longer  the  farther  the  amplifier  is  removed 
from  the  oscillator.  An  amplifier  cell  of  fixed  length  1  appears  to 
have  an  overall  gain  which  is  a  function  of  the  separation  d  from  the 
oscillator.  The  parameter  d  is  clearly  not  an  intrinsic  property  of 
the  amplifier. 

A  satisfactory  explanation  of  this  curious  behavior  can  be  given 
in  terms  of  a  filamentary  structure  of  the  beams.  At  the  output  of 
the  Raman  oscillator  the  filaments  of  Stokes  and  laser  radiation 
coincide  in  space  and  time.  Since  the  Stokes  radiation  has  an  aper¬ 
ture  of  about  10“2  rad,  an  order  of  magnitude  larger  than  the  laser 
beam  aperture,  this  match  between  the  filaments  at  u>l  and  ws  is  at 
least  partially  lost  after  a  distance  d  ~  100a,  where  a  ~  10'2  -  10'1 
cm  is  a  reasonable  dimension  for  the  filament.  Smalle^diameten 
of  the  filament  are  precluded  by  the  observed  angular  definition  of 
Stokes  radiation.  The  gain  of  the  amplifier  depends  on  the  spatial 
distribution  of  the  Stokes  input  with  respect  to  the  filamentary  struc¬ 
ture  of  laser.  One  has  a  typical  insertion  loss  for  the  Stokes  input 
into  each  filament.  The  low  initial  slope  should  correspond  to  an 
average  gain,  as  calculated  from  Eq.  (5),  if  one  determines  |E[J2 
from  the  average  laser  flux  density.  For  benzene  the  low  initial 
amplifier  gain  for  laser  and  Stokes  light  polarized  parallel  to  each 
other  was  found  to  be  gs  ~  0.025  cm  1  for  a  laser  flux  density  of 
107  watts/cm2.  According  to  Eq.  (5),  this  corresponds  to  a  Raman 
susceptibility  of  about  1.2  x  10* 12  esu,  in  good  agreement  with  the 
value  determined  from  Eq.  (13)  for  a  total  Raman  scattering  cross 
section  or  =  5.6  x  10"2S  cm2*  as  determined  by  Damen,  Leite,  and 
Porto,20  and  a  vibrational  width  Tv  =  3.1  cm-1. 

The  gain  in  the  high-intensity  filaments  can  be  higher  by  as  much 
as  a  factor  of  60.  This  enhanced  gain  becomes  evident  when  the 
Stokes  intensity  in  these  channels  has  built  up  to  a  dominant  level. 
The  curves  in  Fig.  2  show  that  the  insertion  loss  into  the  high-gain 
channel  is  about  a  factor  of  10‘2  for  the  geometry  used.  When  the 
mode  structure  of  the  laser  is  cleaned  up,  the  observed  gain  en- 


r 


150 


PHYSICS  OF  QUANTUM  ELECTRONICS 


Banzan*  concentration,  % 

Fig.  4.  The  power  gain  constant  per  unit  length  gs  for  the  benzene 
Stokes  line  as  a  function  of  relative  concentration  in  benzene-nitro¬ 
benzene  (A)  and  benzene-HCe'one  (o)  mixtures. 

hancement  is  reduced  in  agreemenl  with  the  proposed  model. 

The  linear  parts  of  amplifier-response  curves  of  the  type  shown 
in  Fig.  3  may  now  be  used  io  measure  the  Raman  gain  as  a  function 
of  other  variables.  In  Fig.  4  the  gain  per  unit  length  of  benzene 
Stokes  radiation  at  992  cm"1  is  plotted  as  a  function  of  the  relative 
benzene  concentration,  when  the  amplifier  cell  is  filled  with  benzene- 
acetone  and  benzene -nitrobenzene  mixtures.  For  benzene -acetone 
the  gain  is  a  linear  function  of  the  relative  benzene  concentration. 
This  behavior  could  be  expected  if  there  is  little  change  in  damping 
or  local  field  due  to  interaction  of  the  symmetric  vibration  with 
neighboring  molecules.  This  linear  dependence  has  also  been  ob¬ 
served  for  the  spontaneous  emission  of  the  992-cm"1  benzene  line  in 
mixed  liquids.21  This  behavior  should  be  contrasted  with  the  Raman 
oscillator  output  of  mixtures,  where  the  interpretation  is  hazardous 
because  of  nonlinear  competition  between  different  oscillating 
modes.22  The  curve  for  the  nitrobenzene-benzene  mixture  shows  a 
remarkable  behavior.  The  gain  does  not  decrease  down  to  50  per 
cent  benzene  concentration.  The  decrease  in  the  density  of  benzene 
molecules  is  apparently  compensated  by  an  increase  in  the  Raman 
polarizability  per  molecule.  This  may  be  caused  by  stronger  local 
field  corrections  or  by  a  decrease  in  damping.  The  benzene  Stokes 
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gain  is  not  zero  in  pure  nitrobenzene.  This  small  gain  can  be  ex¬ 
plained  by  the  presence  of  a  nitrobenzene  vibration  at  1,000  cm'1. 

It  is  clear  that  the  Raman  gain  can  in  principle  be  measured  for  vi¬ 
brations  that  normally  do  not  show  up  in  oscillators,  because  they 
do  not  have  the  opportunity  to  get  above  threshold. 

Another  example  is  the  measurement  of  the  Raman  gain  as  a 
function  of  polarization.  If  one  inserts  before  the  amplifier  a  quartz 
plate  of  suitable  thickness  and  orientation,  so  that  it  gives  1/4A  re¬ 
tardation  at  oil  and  l/2\  retardation  at  ua,  the  Stokes  gain  for  a 
linear  input  signal  in  the  presence  of  a  circular  polarized  pump  can 
be  measured.  The  input  signal  must  be  decomposed  into  two  circu¬ 
lar  polarized  components.  Since  these  have  different  gain  constants, 
Ksame  and  g0pp’  respectively,  the  output  will  become  elliptical ly 
polarized.  The  polarization  ellipse  of  the  output  can  be  analyzed  in 
a  well-known  manner.  The  ratio  of  the  major  to  minor  axis  is 


exp  <Ksame  0  +  exp  <Sopp  *) 
exP  tesame  f>  ~  exP  <Sopp  f) 


The  experimental  results  give  for  the  ratio  of  the  gain  constants 

-^E.  ,  fb'\  ,  0.44  *  0.05 
gsame  «aT  +  l/! 


A  rotation  of  the  axes  of  the  ellipsoid  would  be  a  measure  of  the 
difference  in  the  real  parts  of  the  Raman  susceptibility  for  the  two 
circular  components.  Although  such  an  effect  has  been  seen,  our 
present  data  are  not  conclusive.  The  real  part  depends  sensitively 
on  the  offset  fiom  exact  vibrational  resonance  of  the  Stokes  compo¬ 
nents. 

When  a  quartz  plate  is  used,  which  gives  A  retardation  at  and 
A/2  retardation  at  wl«  the  linear  polarized  Stokes  input  may  be 
given  an  arbitrary  angle  with  respect  to  the  linear  polarized  laser 
beam  by  rotating  the  quartz  plate.  In  this  manner  it  was  possible  to 
measure  the  ratio  of  parallel  to  perpendicular  gain 


gi 

Ell 


_ 3b'2 

45a'2  +  4b'2 


=  0.19  i  0.07 


The  inversion  and  anisotropy  ratios  are  mutually  consistent.  I  hey 
are,  however,  not  compatible  with  the  depolarization  ratio  deter¬ 
mined  from  the  angular  dependence  of  spontaneous  Raman  scattering 
from  a  polarized  beam  of  a  gas  laser.  Daman,  Leite,  and  Porto20 
give  for  benzene  6b/2/(45a'2  +  7b'2 )  <  0.07.  The  discussion  above 
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under  Polarization  Modes  shows  unequivocally  that  the  same  con¬ 
stants  a'  end  b'  determine  the  polarization  p-operties  of  both  spon¬ 
taneous  and  stimulated  emission.  The  reason  for  this  discrepancy 
is  stiL  obscure,  but  it  will  be  discussed  somewhat  further  in  the 
concluding  section. 


DISCUSSION 

The  material  properties  of  interest  in  the  Raman  effect,  i.e.,  the 
components  of  the  Raman  susceptibility  tensor,  can  be  derived  from 
Raman  amplifier  data.  It  is  possible,  in  principle,  to  measure  sepa¬ 
rately  ihe  effects  of  the  real  and  imaginary  parts.  The  imaginary  part 
can  be  obtained,  probably  more  simply  and  accurately,  from  the 
characteristics  of  spontaneous  Raman  scattering.  The  experiments 
of  Terhune23  on  parametric  generation  of  the  anti -Stokes  frequency 
away  from  the  phase-matched  direction  give  similar  information 
about  the  absolute  value  of  the  Raman  susceptibility. 

The  main  interest  of  the  Raman  amplifier  data,  which  are  gen¬ 
erally  more  useful  for  interpretation  than  oscillator  data,  lies  per¬ 
haps  in  their  information  about  the  properties  of  the  light  beams. 

The  high-intensity  laser  and  Stokes  beams  cann o*  be  considered  as 
homogeneous  distributions  of  constant  intensity.  Even  if  a  laser 
initially  emits  a  monochromatic  diffraction-limited  plane  wave,  such 
a  uniform  structure  may  not  be  preserved  in  the  traversal  of  Bril- 
louin  or  Raman  media.  Chiao,  Garmire,  and  Townes24  have  pointed 
out  the  intensity-dependent  focusing  effect  of  light  beams.  These  ef¬ 
fects  may  be  considerably  enhanced  if  the  real  part  of  the  Raman  or 
Brillouin  susceptibility  near  resonance  is  taken  into  account.  Figure 
3  shows  tiat  the  Raman  gain  in  a  nominal  10  MW/cm2  beam  can  be  a 
factor  e  per  centimeter,  corresponding  to  an  actual  flux  density  in 
the  filament  of  400  MW/cm2.  If  a  Stokes  component  is  1  cm"1  re¬ 
moved  from  exact  resonance  to  the  side  where  Xg  >  the  real  part 
of  the  Stokes  index  of  refraction,  27TXg|ELl2»  proportional  to  the  in¬ 
tensity  of  the  laser  beam,  may  be  as  high  as  10-4.  Since  the  laser 
intensity  may  'ary  appreciably  over  a  lateral  displacement  a  =  1  mm, 
gradients  in  the  index  of  refraction  for  Stokes  light  as  high  as  dn/dx 
«  KT3  cm'1  would  occur.  The  Stokes  light  from  an  area  with  radius 
a  would  be  concentral  d  in  a  distance  f'  ~  (na)^2(dn/dx)'l/2,  or  in  a 
few  centimeters,  if  a  =  1  mm.  This  mechanism  would  favor  the  am¬ 
plification  of  Stokes  light,  displaced  somewhat  to  higher  frequencies. 
As  the  Stokes  intensity  grows  ir.  the  filament,  its  gradient  in  inten¬ 
sity  will  in  turn  focus  the  laser  beam.  It  is  not  clear  that  a  stable 
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steady  situation  with  trapped  beams  will  result.  It  seems  more 
likely,  in  analogy  to  hydrodynamic  instabilities,  that  a  turbulent 
field  is  set  up  by  this  lateral  instability.  The  possibility  exists  that 
the  state  of  polarization  of  the  light  in  this  turbulent  region  is 
changed.  This  might  account  for  the  discrepancy  between  the  ani¬ 
sotropy  of  stimulated  Human  gain  and  the  depolarization  of  sponta¬ 
neous  radiation. 

Turbulent  instability  could  also  account  for  the  stochastic  nature 
of  the  gain  measurements,  for  various  other  characteristics  of  the 
Raman  amplifier  discussed  above,  and  for  Raman  gain  anomalies 
when  a  single-mode  laser  beam  is  shone  into  a  Raman  cell  as  de¬ 
scribed  by  McClung.25  Filamentary  structure  or  multimode  effects 
are,  of  course,  other  terms  to  describe  inhomogencitics  of  the  inten¬ 
sity  in  space  and  time.  At  the  present  time,  the  existence  and  impor¬ 
tance  of  these  nonuniform  intensity  distributions  for  stimulated  gain 
processes  seem  to  be  well  established,  although  their  origin  and 
detailed  characteristics  need  further  clarification. 
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Anomalies  in  the  gain  of  stimulated  Raman 
processes  in  liquids  have  been  investigated 
by  several  groups. 1-3  The  study  with  a  Stokes 
amplifier  cell4  revealed  a  regime  where  the 
gain  per  unit  length  was  one  to  two  orders  of 
magnitude  larger  than  the  value  calculated  from 
the  spontaneous  Raman-emissicn  cross  section. 
The  multimode  theory  of  Bloentbergen  anc;  Shens 
is  inadequate  to  explain  this  discrepancy,  and 
the  experiment  of  McClung6  has  shown  that  the 
same  anomaly  exists  when  the  input  laser  pow¬ 
er  is  essentially  in  a  single  mode.  The  ampli¬ 
fier  cell  studies4  suggested  that  partially  de¬ 
polarized  filaments  of  high  intensity  are  formed 
as  the  laser  beam  passes  through  the  liquid. 
This  self-focusing  action  of  laser  beams  was 
foreseen  by  Chiao,  Garmire,  and  Townes,7 
Askarjan,®  and  Talar.of.®  It  is  due  to  the  inten¬ 
sity-dependent  index  of  refraction  coupled  with 
transverse  gradients  in  the  initial  intensity 
distribution.  This  focusing  action  should  be 
strongest  in  liquids  with  a  large  quadratic  Kerr 
effect  due  to  anisotropic  polarizabilities  of  the 
molecules.®  The  purpose  of  this  Letter  is  to 
describe  new  experimental  results  which  sup¬ 
port  this  view  and  extend  the  considerations 
of  our  very  brief  previous  communication.10 
In  the  meantime  Hauchecorne  and  Mayer"  have 
independently  arrived  at  similar  conclusions  on 
the  basis  of  an  elegant  experiment.  Shen  and 
Shaham1*  have  also  carried  out  experiments 
with  results  similar  to  ours 
In  a  conventional  Raman  oscillator  experi¬ 
ment  with  a  collimated  laser  beam,  one  can 
define  a  threshold  condition  depending  on  pump 


a 
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FIG.  1.  (a)  The  threshold  length  for  stimulated  Stokes 
production  in  a  nitrobenzene  cel!  as  a.  function  of  the 
length  of  a  cell  filled  with  bromobenzene  placed  imme¬ 
diately  in  front.  The  vertical  dashed  line  indicates  the 
threshold  for  Stokes  production  in  bromobenzene, 

(b)  The  threshold  length  for  stimulated  Stokes  produc¬ 
tion  in  a  nitrobenzene  cell  preceded  by  a  30 -cm  long 
bromobenzene  cell,  as  a  function  of  the  distance  be¬ 
tween  the  two  ceils. 
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intensity  and  cell  length,  such  that  i*th^tliv  =  con" 
stant,  where  the  exponent  has  the  experimen¬ 
tal  range  of  value  1  <.v  <  2.  The  threshold  con¬ 
stant  can  be  lowered  significantly,  if  the  laser 
beam  first  traverses  a  cell  with  another  fluid. 
Figure  1(a)  shows  the  reduction  in  the  thresh¬ 
old  length  of  a  nitrobenzene  cell  as  the  length 
of  a  cell  filled  with  bromobenzene  and  placed 
immediately  in  front  of  the  nitrobenzene  cell 
is  increased.  Note  that  the  reduction  occurs 
even  before  the  first  cell  reaches  its  Stokes 
limit.  The  Stokes  shift  of  the  substance  in  the 
first  cell  is  far  from  the  nitrobenzene  Stokes 
shift  of  1345  cm”1.  In  addition,  filters  are 
used  to  insure  complete  decoupling  between 
the  cells.  The  Stokes  emission  in  the  first 
cell  may  also  be  suppressed  by  selective  ab¬ 
sorbers  without  affecting  the  creation  of  hot 
filaments.  The  reduction  of  a  factor  10  in  thresh¬ 
old  length  means  that  the  effective  pump  inten¬ 
sity  is  enhanced  by  almost  two  orders  of  mag¬ 
nitude.  If  the  distance  between  the  two  cells 
is  increased,  the  reduction  is  decreased.  The 
effect  of  the  first  cell  does  not  extend  beyond 
20  cm  from  its  end  face,  as  shown  in  Fig.  1(b). 

We  have  also  shown  the  enhancement  of  inten¬ 
sity  in  filaments  by  means  of  another  nonlinear 
process,  the  harmonic  production  in  a  quartz 
platelet  placed  at  the  end  of  the  cell.  This  tech¬ 
nique  has  been  used  in  a  more  elegant  and  quan¬ 
titative  manner  by  Haucheeorne  and  Mayer.11 

The  creation  of  the  filaments  in  the  first  cell 


has  been  recorded  on  a  Polaroid  film.  The 
difficulty  in  making  them  visible  for  the  cell 
lengths  used  iri  our  experiments  is  that  the  main 
fraction  of  the  laser  intensity  goes  through 
without  taking  part  in  the  focusing  process. 

Since  we  had  previous  evidence  that  the  effec¬ 
tive  pump  light  is  partially  depolarized,  we 
put  the  cell  between  a  crossed  polarizer  and 
analyzer  to  suppress  most  of  the  laser  light. 

The  image  of  the  end  face  of  the  cell  is  record¬ 
ed  in  Fig.  2.  When  the  cell  length  is  not  suf¬ 
ficient  for  self-focusing,  a  weak  homogeneous 
intensity  distribution  shown  in  Figs.  2(a)  and 
2(b)  appears.  When  self-focusing  occurs  with¬ 
in  the  cell  length,  bright  spots  appear.  Their 
num  ,er  Increases  with  increasing  cell  length. 
Their  diameter  is  between  20  and  80  p.  This 
dimension  is  in  agreement  with  the  loss  of  the 
filamentary  structure  due  to  diffraction  in  20 
cm  of  air.  The  pattern  changes  from  laser 
shot  to  laser  shot,  unless  strong  reproducible 
inhomogeneities  are  introduced  into  the  beam 
by  apertures.  Figure  2(e)  shows  that  filaments 
preferentially  form  in  a  circular  pattern  after 
diffraction  of  the  beam  through  a  circular  hole 
before  entering  the  cell. 

There  is  a  very  good  correlation  between 
the  self-focusing  ability  of  a  liquid  and  its  in¬ 
tensity-dependent  index  of  refraction.9  Liquids 
are  decreasingly  effective  in  the  order  CS2, 
nitrobenzene,  bromobenzene,  benzene,  and 
act  tone,  while  no  detectable  focusing  action 


FIG.  2.  Photographs  of  the  end  face  of  a  nitrobenzene  cell,  between  crossed  polarizer  and  analyzer,  of  variable 
length  indicated  in  mm.  (a)  and  (h)  The  length  is  below  threshold,  (c)  The  iength  is  at  threshold  for  Stokes  produc  - 
tion.  (d)  The  length  is  well  above  threshold,  (el  The  laser  beam  is  diffracted  by  a  small  circular  aperture  before 
entering  the  cell.  Filaments  are  found  on  rings  of  maximum  transverse  intensity  gradient. 
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occurred  with  isopentane,  alcohol,  water,  and 
CCl4.  We  have  found  no  good  correlation  with 
the  electrostrictive  properties  nor  with  a  reso¬ 
nance  associated  with  two-quanta  absorption. 
Theory  also  predicts  that  their  contribution 
to  the  intensity-dependent  index  should  be  small¬ 
er  than  that  of  the  anisotropic  polarizability. 

The  experimental  threshold  for  Raman  laser 
action  appears  to  bo  determined  not  so  much 
by  the  value  of  the  Raman  susceptibility,  but 
rather  by  the  self-focusing  capability  of  the 
fluid.  This  suggests  that  stimulated  Raman 
oscillations  may  be  induced  in  liquids  that  re¬ 
main  below  threshold  in  the  unfocused  laser 
beam  by  mixing  them  with  a  focusing  liquid. 

A  cell  with  pure  cyclohexane  was  below  thresh¬ 
old  in  the  collimated  laser  beam.  Stimulated 
Stokes  lines  of  cyclohexane  were  produced  if 
the  cell  was  filled  with  a  mixture  of  two-thirds 
cyclohexane  and  one-third  nitrobenzene.  The 
totally  symmetric  CC14  molecule  does  not  have 
sufficient  focusing  action  of  its  own.  We  have 
observed  three  orders  of  stimulated  Stokes 
light  with  the  characteristic  shift  of  459  cm-1 
of  the  totally  symmetric  vibration  of  CC14 
in  a  mixture  of  75  cc  CC14,  15  cc  nitroben¬ 
zene,  and  10  cc  CS2.  It  was  necessary  to  use 
a  mixture  of  the  last  two  focusing  fluids  in  or¬ 
der  to  stay  below  the  Stokes  threshold  of  both 
of  these. 

The  characteristic  focusing  length  follows 
from  the  curvature  of  light  rays  in  a  transverse 
gradient  of  index  of  refraction.  This  length 
should  be  smaller  than  the  Fresnel  length  of 
the  beam,  d2/\,  and  than  the  cell  length.  If 
the  intensity  drops  by  an  amount  S0  over  a  dis¬ 
tance  a  on  either  side  of  an  intensity  maximum 
in  the  beam,  the  self-focusing  length  is 


where  m20o  is  the  intensity-dependent  change 
in  index  of  refraction.  A  more  detailed  calcu¬ 
lation  of  the  focusing  action  in  the  laser  beam 
has  recently  been  given  by  Kelley.13  If  one 
takes  a  equal  to  the  size  of  the  laser  beam  d 
and  puts  g0  equal  to  the  Intensity  at  the  beam 


center,  he  finds  that  the  focusing  length  under 
the  conditions  of  our  experiments  is  1  m  for 
CS2  with  m2  10"m  esu;  and  3  m  for  benzene 
with  m2  =  0.2  x  10-n  esu  The  experimental 
focusing  lengths  are  an  order  of  magnitude 
smaller.  Furthermore,  the  whole  beam  does 
not  focus  as  a  unit.  The  multimode  effect  must 
still  be  invoked  to  explain  this  discrepancy. 
Apparently  random  small-scale  intensity  varia¬ 
tions  present  in  the  input  laser  signal  are  re¬ 
sponsible  for  more  rapid  and  piecemeal  focus¬ 
ing  of  the  beam.  It  is  not  clear  at  this  moment 
what  factors  determine  the  diameter  of  the 
filaments.  Apparently  some  opposing  defocus- 
ing  effects  set  in  after  the  filaments  have  reached 
a  given  size  or  intensity. 
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